Graphene nanoribbons subject to gentle bends 
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Since graphene nanoribbons are thin and flimsy, they need support. Support gives firm ground 
for applications, and adhesion holds ribbons flat, although not necessarily straight: ribbons with 
high aspect ratio are prone to bend. The effects of bending on ribbons' electronic properties, 
however, are unknown. Therefore, this article examines the electromechanics of planar and gently 
bent graphene nanoribbons. Simulations with density-functional tight-binding and revised periodic 
boundary conditions show that gentle bends in armchair ribbons can cause significant widening or 
narrowing of energy gaps. Moreover, in zigzag ribbons sizeable energy gaps can be opened due 
to axial symmetry breaking, even without magnetism. These results infer that, in the electronic 
measurements of supported ribbons, such bends must be heeded. 
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I. INTRODUCTION 

Graphene nanoribbons (GNRs) are atomically thin 
and only nanometers wide, which makes them the flim- 
siest materials in the world. Today such ribbons, ac- 
claimed fo r p romising applications, are fabricated in 
many ways^l. and investigated for heat conduction 3 , 
edge features-^H and electronic characteristics-^, among 
many other properties. However, since the ribbons are 
flimsy, they need stabilizing support — although even then 
ribbons can get folded, torn, rippled, and bentP^^ 

Also supports are different, as the interaction with 
graphene can be either physical or chemical. In physi- 
sorption the support interaction is weak, graphene's elec- 
tronic structure remains decoupled, and adhesion arises 
from the dispersive van der Waals interactions aloneP^ 
In chemisorption the support interaction is stronger, and 
the presence of chemical bonds alters graphene's elec- 
tronic structured Therefore adhesion, responsible for 
holding ribbo ns pla nar, ranges from e = 4 meV/atom to 
70 meV/atomP^L^ However, fabrication processes, sur- 
face inhomogeneities, pinning, AFM tip manipulation, 
heat treatment, or mechanical strains can make ribbons 
subject to gentle bends, as sketched in Fig. ITJa). In- 
deed, planar and gentle bending can be directly seen in 
scanning tunneling microscopy experimentsPHSS Distor- 
tions like twisting, on the contrary, are less relevant on 
supportsP^EIDOnly gentle bends are interesting, as sharp 
bends are struct urall y unstable (ribbons would desorb 
and fold instead) PS] 

The purpose of this work, therefore, is to answer the 
following simple question: What happens to GNRs' elec- 
tronic structure upon planar bending? It turns out 
that simple geometrical arguments, together with nearest 
(and next-nearest) neighbor tight-binding reasoning, are 
sufficient for a thorough understanding of the electrome- 
chanics of bent GNRs. These insights should hence help 
interpreting imperfect experiments with these distortion- 
prone ribbons. 



II. SIMULATING PHYSISORBED RIBBONS 
WITH PURE BENDING 

I modeled GNRs as free-standing, without explicit 
presence of the support; it was there merely as a planar 



constraint. The underlying justification was to model ph- 
ysisorption where the support and GNR electronic struc- 
tures are essentially decoupled. For chemisorption the 
results are not directly valid. 

The focus is on the bent sections of very long rib- 
bons, with bending viewed as a local property. Apart 
from bending, planar ribbons can also stretch and 
shear; those deformations have been investigated by 
conventional methods ! 17 * 18 ! The deformation mode cer- 
tainly depends on the experimental conditions, and es- 
pecially in s hort r ibbons the strain patterns can become 
complicatedf^Ell However, ribbons yield easily upon lat- 
eral forcing, and adjust themselves readily to minimum- 
energy geometries J22H221 l n bent geometries sliding is par- 
ticularly easy as the ribbon and the support are mostly 
out of registry. Long ribbons pinned at two distant lo- 
cations, therefore, can remove high-energy shearing and 
stretching by sliding, and favor pure bending! 2 ^ I remark 
that the central results, as discussed below, will be valid 
also beyond pure bending. 

I modeled the electronic structure by density- 
functional tight-binding (DFTB) method 221221 , and the 
bent geometry itself by revised periodic boundary 
conditions. 28 " 30 -^ Atoms in the simulation cell were from 




FIG. 1: (Color online) (a) Supported graphene nanoribbon 
sketched for gentle, planar bends, (b) The red (dark gray) 
atoms constitute the unit cell, simulated with revised peri- 
odic boundary conditions; the symmetry operation is rota- 
tion of an angle a around given origin (not shown). Ribbon's 
width W is defined by the outmost carbon atoms and R is the 
mean radius of curvature. The bending parameter of Eq.(|TJ, 
here Q = 0.1, also equals (approximately) to the compressive 
strain at the inner edge (ei n = —6) and to the tensile strain 
at the outer edge (e ut = 6). 
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Bending parameters 

FIG. 2: (Color online) The cross-ribbon averaged strain e a v g 
as a function of bending parameter for AGNRs (solid lines) 
and for ZGNRs (dashed lines); the bold dashed line is an an- 
alytical estimate. Inset: Elastic energy density as a function 
of bending; the bold dashed line is an analytical estimate, 
Eq.Q. In both plots the line width is proportional to W. 
Hence both in e avg and in E/A the largest deviations are for 
the narrowest ribbons. 



III. BENT RIBBONS GET STRETCHED 

Let us now turn attention to the results. Prior to dis- 
cussing electronic properties, however, let us first look 
at energy and geometry. The energy in bent ribbon, as 
hinted by nanoshell elasticity^, comes chiefly from axial 
in-plain strain. A quick estimate yields energy per unit 
area as E/A = gfc9 2 , where k = 25 eVA -2 is graphene's 
in-plane modulus.^ This simple estimate is in fair agree- 
ment with the simulations, as shown in the inset of Fig. [2] 
Only the narrowest ribbons deviate from this estimate, 
for two reasons: First, the comparison of W between 
atomic and continuum methods is inherently ambiguous; 
for small W this ambiguity is emphasized. Second, in 
narrow ribbons the value of k is affected by distinct elas- 
tic properties near the edges. While I could remedy these 
deficiencies by improving the model, my main interest is 
not in the minutiae of narrow ribbons, but in the wider 
ribbons and their universal trends. 

If we set the adhesion energy e equal to the strain en- 
ergy, we get 



GNR translational cell of length L, and the associated 
symmetry operation was a rotation of an angle a around 
a given origin, as in Refs. I2U1 and I3T1 [see Fig.JTJb)]. This 
means, therefore, that the simulated systems were effec- 
tively GNR hoops, containing hundreds of thousands of 
atoms. Since the bends are gentle and the charge trans- 
fer with physisorption usually small, it's reasonable to 
assume that simulation describes the properties of bends 
in GNRs in a local senseP At any rate, regarding the 
bending, simulations were exact and the sole approxima- 
tion was the DFTB method itself. 

Throughout this article I will use the dimensionless 
parameter 



9 £ = y/6e/kA c 
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to quantify the amount of bendingf^I Then, to simu- 
late GNR of width W with bending close to 9', I chose 
a = L/R', with R' = W/(2&) as the initial guess for the 
radius of curvature, and optimized the structure. The 
only fixed parameter was a, so R and 9 were outcomes 
of the optimization, although R sa R 1 and 6«0'. Here 
I remark that, because a's are small (down to 1CP 3 ra- 
dians), the optimization was arduous and required max- 
imum force criteria as small as / lnax < 10 _ft eV/A (look 
Ref. [30] to see why) . 

I conducted such simulations for hydrogen-passivated 
armchair ribbons (7V-AGNRs) and zigzag ribbons (N- 
ZGNRs), with N = 5 . . .40, with W up to 84 A, with 74 
different GNRs in total (see Ref. [33J for GNR notations) . 
Each GNR was optimized for 10 bendings between = 
and 8 = 0.1. The reason for 9 = 0.1 as the upper limit 
for bending that I term "gentle" will be clarified later. 
Finally, the number of K-points was 50 k/L for geom- 
etry optimization and 500 k/L for electronic structure 
analysis. 



as a rough estimate for the limit where the ribbon 
rather desorbs and straightens than remains adsorbed 
and bent on the support. The maximum adhesion e = 
70 meV/atom yields 9 C = 0.08, justifying the upper limit 
9 < 0.1 for a "gently" bent ribbon, even though 9 e re- 
ally depends on the substrate. This is only an order- 
of-magnitude estimate, as fluctuations and finite-size ef- 
fects can cause ribbons to desorb earlier. Direct experi- 
mental evidence^ shows how GNRs on Si02 bend up to 
9 sa 0.01 — still nearly half the simple-minded limit of 
9 e ~ 0.025 given by e « 6 meV/A 2 P 

Given the definition for 9, the strain on ribbon's in- 
ner edge is e- m « —9 and on the outer edge e out ~ 9. 
With strains around 10 % the bond anharmonicities be- 
gin to emerge, and stretching becomes cheaper, compres- 
sion more expensive. This implies that the neutral line 
moves away from the origin (i? increases) and ribbon 
stretches. We can take this effect into account by a 
strain-dependent in-plane modulus, k{e) = fco(l — 75). 
Then, by minimizing the total energy per unit length 



R+W/2 



R-W/2 



-k (l-je)e 2 dr 



(3) 



with respect to R, we obtain the cross-ribbon averaged 
strain as 



;7@ 



(4) 



This analytical estimate, given the value 7 = 1.7 ob- 
tained from DFTB simulations of stretched GNRs, agrees 
well with simulations, as shown in Fig.j2l The largest de- 
viations occur again for the narrow ribbons, albeit with 
opposite tendencies for AGNRs and ZGNRs due to dif- 
ferent edge morphologies. 
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FIG. 3: (Color online) Electronic structure of bent AG- 
NRs. (a) Bending-induced gap changes E g (Q) — E g (0) for 
the three q- families of AGNRs. Line width is proportional 
to W; dashed lines are estimates for q = 0, 1. Inset: Gaps 
in straight AGNRs. (b) Wave functions for the frontier or- 
bitals in straight 16-AGNR (g = 1): the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied molec- 
ular orbital (LUMO). (c) Density of states for 16-AGNR 
with straight (O = 0.0), bent (O = 0.1), and stretched 
(s = §7(0.1) 2 = 0.85 %) geometries. 



along the ribbon's axis, and therefore stretching tends to 
narrow the gap (for q = AGNRs the situation is the 
opposite and for q = 2 intermediate).^ Next, if we pre- 
tend, in effect, that the bent AGNRs experience only the 
average axial strain (even if the strain is uneven), and 
thus juxtapose e with e avg from Eq.(|4l, we get 



A£ 9 (e) « -(-i)ve 2 



(5) 



Figure [3la) plots these estimates for q = and q = 1 
AGNRsby the dashed lines. The fair agreement sug- 
gests that the electronic structure of AGNRs subject to 
bending is dominated by the cross-ribbon averaged strain. 
Similar physics has been observed p reviously i n bent car- 
bon nanotubes and twisted GNRsP^ 31 ' 39 ' 40 ' Hence the 
argument is easily generalized to combined bending and 
stretching, where the electronic structure is modified by 
the average strain e s trctch + ^jQ 2 ^ 

These trends, as given by four-valence DFTB, are re- 
produced by a 7r-only tight-binding Hamiltonian 



H = 



next n.n 



eje 



j • 



(6) 



with the nearest-neighbor (n.n.) hopping parameter 
t(r) = 2.6 eV - 5.8 eV/A(r - 1.42 A) 



(7) 



and with the next-nearest neighbor hopping equal to zero 
(<' = 0, not shown). Figure [3lc) shows further that 
the stretching analogy extends Deyond energy gaps, as 
the entire density of states (DOS) is well described by 
the stretched geometry. This carries the average strain 
analogy also for optical transitions, as shown earlier.^ 
Note that, if gauged through the relative gap change 

\AE g /E g \ w 0.8k~ 1 W@ 2 , the influence of bending be- 
comes more important as W increases. 



IV. ARMCHAIR RIBBONS ARE DOMINATED 
BY STRETCHING 

Equipped with these geometrical notions, let us now 
turn attention to the electronic properties, starting with 
AGNRs. The inset in Fig. [3]ja) shows the energy gaps 
for the known three families of A^-AGNRs, defined by 
q = mod (AT ; 3). I3313U The gaps scale as Eg ~ {3W~\ 

where /? rj 13 eVA for q — 0, 1 (for q = 2 the scaling 
is bit different). Figure pFa) shows how these gaps re- 
spond to bending: they widen or narrow with the same 
y-dependent families. Deviations occur only for the nar- 
rowest ribbons. 

These trends can be understood by the following 
model. The energy gaps in stretched q = 0, 1 AG- 
NRs depend on the strain as A^ trai s ht rj (-1)*e6 with 
5 = 12 eV (fit for q = 2 is just more complex) ! 17 l 38 l The 
origin for this strain-dependence is illustrated in Fig|3[b) 
for 16-AGNR with q = 1: the highest occupied orbital is 
bonding and the lowest unoccupied orbital is antibonding 



V. ZIGZAG RIBBONS ARE DOMINATED BY 
BROKEN SYMMETRY 

Let us now leave AGNRs and turn our attention to the 
electronic properties of ZGNRs. First I have to remind 
that, for straight ZGNRs as such, the spin-parallel DFTB 
simulations are dubious, given the prediction for a spin- 
polarized ground state! 3 - 3 -' The magnetic structure should 
arise when the curious flat bands near the Fermi-leveP^ 
[left panel of FigQa)] — the famous edge states — spin- 
polarize, lift degeneracies, and open a gap (not shown) 
This way spontaneous magnetization can stabilize the 
electronic structure. 

It has been shown that, unlike in AGNRs^, 
the electr onic structure in ZGNRs is unaffected by 
stretching! 17 ' 42 ' Therefore, after discovering above the 
average-strain argument with AGNRs, it's natural to 
guess that bending would leave ZGNRs' electronic struc- 
ture unaffected. The right panel of Fig.Qa) shows, how- 
ever, that bending can open an energy gap in ZGNRs. 
That is, electronic structure is stabilized by sheer bend- 
ing, and the cause for magnetic spin-polarization is lost. 
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The Hamiltonian ([6]), with hoppings Q and Q, re- 
produces the electronic structure fairly well, as shown 
by the DOS for 10-ZGNR in Fig. Sjc). [Pure stretch- 
ing leaves DOS intact (not shownjj. Third step; the 
flat band energies (also the band dispersion) are propor- 
tional to t' and hence proportional to edge strain via 
Eq.j8|. Upon bending, the reflection symmetry breaks 
and states localize on either of the edges with strain dif- 
ference e ut — £in = 28; opposite edges hence get unequal 
hoppings At' — t' out — t' in oc O. Because energy splitting 
is proportional to At', it is also proportional to O. This 
is the mechanism how bending splits the flat bands with 
direct proportionality to O. 

As mentioned above, since t' gives flat bands a small 
dispersion, splitting does not open the gap immediately. 
When W increases, the span of the flat region in fc z -space 
increases, and gap opening requires larger splitting. A fit 
to all ZGNRs yields a critical value for opening a gap as 
e crit « W/200 nm (or i? crit 100 nm for all W), yielding 
the energy gap as 



£ ZGNR ^ 4 eV( Q _ @ cr . t) (Q) 



The gaps, displaying values up to 0.4 eV, are plotted in 
the inset of Fig. Etc). 



FIG. 4: (Color online) Electronic structure of bent ZGNRs. 
(a) Band structure of 10-ZNGR with a straight and bent ge- 
ometries. For the straight ribbon, having a reflection symme- 
try, the dashed lines denote symmetric and solid lines denote 
antisymmetric states under reflection. (For the bent ribbons 
no such distinction can be made.) (b) Wave functions of 
the frontier orbitals in straight 10-ZGNR: highest occupied 
molecular orbital (HOMO, symmetric) and lowest unoccu- 
pied molecular orbital (LUMO, antisymmetric), (c) Density 
of states in straight and bent 10-ZGNR. The dashed lines 
are from a next-nearest neighbor tight-binding model [Eq. Q 
with Eq.jljl]. Inset: Gaps for all ZGNRs plotted as a function 
of 9 - Ocrit (Eg = when 9 < G cr it). 



Note how the band structure changes only near the 
Fermi-level, while other bands remain stable. 

The mechanism of the gap opening is related to broken 
reflection symmetry, as clarified by the following three- 
step reasoning: First step; the bands are called "flat" 
because they have small dispersion. In the absence of 
next-nearest neighbor hopping (t' = 0), the Hamilto- 
nian ^ gives edge states whose dispersion and energy 
are essentially zero, independent of t. This is illustrated 
in Fig. pfb), where flat band electrons appear localized 
to next-nearest neighbor sites, separated by vacancies. 
Therefore the t in Eq. Q, even if strain-dependent, 
doesn't affect the flat bands — splitting and dispersion 
hence require next-nearest neighbor hopping t' . Second 
step; fitting t! to strained GNRs by DFTB gives 

t'(r) = 0.25 eV - 0.6 eVA(r - 2.46 A). (8) 



VI. CONCLUSIONS 

The physics in AGNRs and ZGNRs appear hence 
quite different: AGNRs are governed by average strain, 
whereas ZGNRs are governed by broken reflection sym- 
metry. The effects of broken symmetry on AGNRs or 
average strain on ZGNRs surely exist, but they are just 
less important. In ZGNRs bending can have particular 
impact on transport, since the localization of edge states 
depends on the direction of bending; if ribbon has bends 
both to the left and to the right, the current-carrying 
electrons need to jump from one edge to the other, sug- 
gesting width-dependent resistivitypSl Although it's plau- 
sible that bent ZGNRs indeed acquire gaps and turn 
nonmagnetic, spin-polarized calculations would be op- 
portune, even if the existence of magnetism has been 
disputed also for the straight ribbons 

I obtained similar results also for unpassivated GNRs, 
observing similar phenomena. Thus it appears that these 
clear trends arise from simple physics with plausible ex- 
planations, and it's unlikely that, say, higher level elec- 
tronic structure methods should change the picture. I 
believe, therefore, that these general trends are helpful 
enough to serve as rules of thumb to aid GNR device 
fabrication and analysis. 
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